1 

Stable XOR-based Policies for the Broadcast 
Erasure Channel with Feedback 

Sophia Athanasiadou, Marios Gatzianas, Member, IEEE, Leonidas Georgiadis Senior Member, IEEE 

and Leandros Tassiulas, Fellow, IEEE 



(N 



Abstract 

In this paper we describe a network coding scheme for the Broadcast Erasure Channel with multiple stochastic 
flows, in the case of a single source transmitting packets to N users, where per-slot feedback is fed back to the 
transmitter in the form of ACK/NACK messages. This scheme uses only binary (XOR) operations and involves a 
{SJ . network of virtual queues, along with special rules for coding packets and moving them among the queues that 

ensure instantaneous decodability. Based on the stabilizing scheduling policy described in the framework of (TJ, we 
show that, for N = 4 and i.i.d. erasure events, the network stability region of such a system effectively coincides 
" with its information-theoretic capacity region, and provide a stabilizing policy that employs this XOR-based scheme. 

o 

I. Introduction 

The information-theoretic capacity region of the Broadcast Erasure Channel (BEC) in the case of one transmitter 
c/2 . and iV multiple unicast sessions has been recently determined in O and Q. Both of these papers proposed 
i capacity achieving algorithms based on transmission of linear combinations of packets. However, these schemes 

■ are characterized by high complexity (as operations take place in a sufficiently large sized finite field) and decoding 
^ , delay, since a sufficient number of linear combinations has to be received until a packet is decoded. In flU, we 
00 [ proposed a network coding scheme that overcomes these obstacles by using only XOR operations, generalizing the 

■ 2-user network coding scheme in [5] to the case of 3 users. Thus, two low complexity algorithms were proposed, 
^ . namely X0R1 and X0R2, that additionally had the advantageous property of "instantaneous decodability". By this 

• | term, it is meant that a receiver is able to decode packet p destined for it as soon as it receives an XOR combination 

■ of packets containing p. Algorithm X0R2 was proved to achieve capacity for the case of i.i.d. channels as well as 

■ spatially independent channels with erasure probabilities that do not exceed 8/9. 

However, the system considered in (4| is a saturated system, where a predefined number of packets needs to 

■ be transmitted to each user. This model is not frequently encountered in practice. Moreover, algorithms X0R1 and 
^ , X0R2 cannot be easily generalized to more than 3 users. This happens because, at each time slot, coding choices 

$— 1 1 have to be rigorously determined so that each transmission is optimally exploited in terms of allowing multiple 
. 5^ ■ users to simultaneously decode their packets as well as create favorable future coding opportunities. However, for 
N > 3, the number of coding choices increases dramatically so that there is no clear intuition on the optimal 
choice. 

In the current work, we propose a general network coding scheme for the case of a single transmitter sending 
packets to N users through the BEC with feedback, generalizing the scheme proposed in |4). Any packet arriving 
to the transmitter is placed in one of N queues. Depending on the received feedback, these packets (or XOR 
combinations of them) may travel through a network of virtual queues, before they reach their destination, in order 
to exploit the overhearing benefit of the broadcast channel. Coding and packet movement rules are imposed in order 
to ensure instantaneous decodability of packets and better exploitation of coding opportunities. 
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Sciences, EPFL, Switzerland. L. Tassiulas is with the Computer Engineering and Telecommunications Department, University of Thessaly, 
Volos 38 221, Greece. 
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While in |4| we examined a saturated system, in this paper we consider the more realistic model where packets 
may arrive randomly at the transmitter at any time slot. Additionally, we use a backpressure type online algorithm 
that makes each coding choice based on instantaneous quantities, such as queue sizes, without requiring knowledge 
of future events. Therefore, we do not need to predefine the coding choices, and the proposed network coding 
scheme can be applied for an arbitrary number of users. For the specific case of 4 users and i.i.d. erasure events, 
we present a stabilizing policy on top of this network coding scheme and prove that the policy stability region 
coincides with the information theoretic capacity region of the standard BEC with feedback. This result is quite 
intriguing, considering that the proposed policy only uses a subset of all possible coding choices. Moreover, we 
address inplementation issues and present an overhead analysis for the maximum number of packet addresses that 
need to be included in the header of an XOR-coded packet. 

The network stability of single hop broadcast erasure channesl with feedback has also been examined in |6j, 
which considered broadcast traffic only and investigated the stability regions of plain retransmission and linear 
network coding schemes (parameterized over the field size) as opposed to a proposed dynamic virtual queue-based 
policy. The latter policy was shown to be optimal for 2 users while, for N > 2 and i.i.d. erasures, it achieved a 
stable rate that differs from the cut-set bound by a factor of 0(e m+1 ), where m is the number of queue "levels" 
that participate in the coding decision (see for more details and definitions; m can be loosely regarded as a 
measure of the encoding complexity) and e is the erasure probability. Although the structure of the virtual queues 
and coding rules are inspired by similar concepts as in our work, the actual rules for moving packets between the 
queues are much more involved in our work since we are interested in achieving the optimal stability region for all 
values of e instead of only asymptotic optimality as e — > (these notions of optimality ignore any overhead). An 
additional cause for rule complexity in our work is the fact that multiple unicast sessions are much more difficult 
to handle (due to the inherent competition between different sessions) than a single broadcast session. Furthermore, 
there is no guarantee in O, for the general case of N users, regarding instantaneous decodability. 

The work in Q studied a network which is described by an underlying complete graph where each edge is 
modeled as a Markov chain ON/OFF channel (i.e. a generalization of the memoryless erasure channel), while 
there also exists a special "relay" node with XOR coding capabilities which can overhear all transmissions. Any 
transmissions to/from the relay are error-free. The work considers multiple unicast flows, originating in all nodes 
except for the relay, and explicitly accounts for instantaneous decodability by mapping this constraint into a specially 
constructed conflict graph (a similar graph structure is used in [?] to model the same constraint). It proposes an 
online backpressure policy that requires computing in each slot the maximum weight independent set of the time- 
varying conflict graph. Although the work bears similarities to our paper in terms of mathematical techniques and 
the optimization problem that results, the model is quite different. Hence, the proposed coding policies are quite 
different and the results in cannot be used to show one of our main results, namely that the proposed scheduling 
and coding policies achieve channel capacity for BEC with i.i.d. erasures. In particular, the broadcast channel at 
the relay (which is the only node that can perform XOR coding) is error-free in Q, while we are interested in 
broadcast erasure channels. 

In summary, the contribution of this paper lies in the development of a systematic framework for constructing 
instantaneously decodable feedback-based XOR coding schemes for the BEC with multiple unicast traffic and 
an arbitrary number of users. Although this framework provides general guidelines and determines the required 
properties of such a scheme, it does not specify all details of the coding procedure, thus allowing the code designer 
to trade complexity for performance. Additionally, combining this framework with the backpressure stabilizing 
policy in the model of HI, and carefully limiting the allowable coding choices, leads to the unexpected result that 
proper XOR combining is sufficient to achieve the capacity of the 4-user BEC with i.i.d. erasures. 

The rest of the paper is organized as follows: in Section [TTJ the system model is introduced along with some 
useful notation. In Section [Till the proposed network coding scheme is described, while in Section [TV] the applied 
stabilizing policy is presented. In Section [VJ an outer bound on the stability region of the system under study is 
derived. In Section [VTJ we prove, for the case of 4 users and i.i.d. erasure events, that the stability region of such a 
system coincides with the capacity outer bound of the standard broadcast erasure channel with feedback. In Section 
IVIII we examine some implementation issues while Section I VIII I concludes the paper. 
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II. System model and notation 

We describe some notation that will be used in the following. Sets are denoted by calligraphic letters, e.g. 
A4, and the empty set by 0. The cardinality of set M. is denoted by \M\ and we write M = \M \ . Random 
variables are denoted by capital letters and their values by small case letters. Vectors are denoted by bold letters, 
e.g. A = (Ai, . . . , A n ). The expected value of a random vector is the vector consisting of the expected values of 
its components, i.e., E [A] = (E (Ai) , . . . , E (A n )) . 

We consider a time-slotted system where slot t = 0, 1, . . . corresponds to the time interval [t, t + 1). The system 
consists of a base station B and a set N = {1,2, . . . ,N} of receivers (users). At the beginning of slot t, A{(t) 
data packets arrive at B with E(Ai) = A»; these packets must be delivered to receiver i and are referred to as 

"flow i" packets, where we denote A(t) = (Ai(t), . . . , An {£)) , t = 0, 1, All packets consist of L bits with 

an individual transmission-time duration of 1 slot. A packet transmitted by B may be either correctly received or 
completely erased by any receiver (broadcast medium). After each transmission, the receivers send feedback to B 
(through an error-free zero-delay channel) informing whether the transmitted packet has been correctly received or 
not (ACK/NACK feedback). We also assume that if no packet is transmitted in a slot (idle slot), then all receivers 
realize that the slot is idle. 

Packet arrivals are assumed to be independent and identically distributed across time, but arbitrarily correlated 
at a given time. That is, the process {A(i)}^ consists of i.i.d. random vectors, while the components of each 
vector A (t) may be arbitrarily correlated. Similarly, packet erasures are i.i.d across time and are initially assumed 
to be arbitrarily correlated at a given time (we later concentrate on the special case of spatially i.i.d. erasures). The 
packet arrival and erasure processes are independent. For subsets S, Q C J\f with Sr\Q = 0, we denote by Ps,g the 
probability that a transmitted packet is erased at all receivers in S and received by all receivers in Q (no condition 
is imposed on packet reception or erasure for receivers in M — (S n Q) ). We also denote by e$ the probability that 
a transmitted packet is erased by all receivers in S, i.e., eg = Ps,&- For simplicity, we slightly abuse the notation 
and write e« or instead of or jy, respectively. 

III. Network coding scheme description 

Exogenous packets arriving at B and being intended for user i G J\f are called "native packets for i". A packet 
is simply termed "native" if it is a native packet for some user (due to the unicast traffic, a packet is native for 
exactly one user). According to the policies to be described below, all transmitted packets are either native, or XOR 
combinations of native packets. In other words, any transmitted packet p can be written as p = ©™ = i Qh where qi 
are native packets and we say that "p contains q" . It is possible, and actually beneficial, for p to contain native 
packets for more than one user. To shorten the description in the following, we say that a packet p is an XOR 
combination of native packets even when p consists of a single native packet. The following definitions will be 
crucial in the subsequent analysis. 

Definition 1. User i is a listener of a packet p iff both of the following conditions are true: 

1) p is an XOR combination of packets, not necessarily native, that i has correctly received. 

2) p contains no native packet for i that is unknown to i. Equivalently, if p contains a native packet q for user 
i, then the packet q is known to (i.e. has already been decoded by) i. 

Notice that the second condition of Definition Q] does not assert that p always contains a native packet q for user 
i, only that the existence of such a packet implies that q is known to i. To illustrate this definition, let a packet p 
be of the form p = p\®p2 where pi, p2 are native packets (not necessarily intended for user i € J\f). Then, i is a 
listener for p if either p\ and p2 have both been received by i, or if p itself has been received by i and any of p\ 
and p2 that may be native packets for i are already known to i. 

Definition 2. User i is a destination of a packet p if either p is a native packet intended for user i that has not been 
received yet by (i.e. is unknown to) i, or if p can be decomposed as an XOR combination of the form p = q © c 
where 

1) q is a native packet intended for i and unknown to i, and 

2) i is a listener of c. 
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We hereafter use the term destination to exclusively refer to the above technical definition. Notice that the 
decomposition of a packet p with destination i alluded to in Definition [2] is unique, since c cannot contain an 
unknown native packet for i due to the second condition of Definition Q] (recall that i is also a listener of c). We 
also denote the unique q in p as q = p(i) and call p(i) the "unknown packet" of i in p. As will be seen, transmitted 
packets may have several receivers as destinations or listeners. As a final note, it is apparent that the notions of 
destination and istener are time-varying, as each user receives new packets over time. 

The main features of the policies to be presented are as follows. 

Basic Algorithmic Features 

1) Any transmitted packet is an XOR combination of native packets. 

2) If a transmitted packet (equivalently, an XOR combination) p contains a native packet q intended for i G N, 
which is unknown to i, then i is a destination for p. Hence, since according to the definition of destination 
it holds either p = qoxp = q@c, where c is known to i, user i can immediately decode q as q = p © c. 

Under the proposed policies, packets may be placed in various virtual queues, based on the received feedback. A 
general virtual queue Qffy is characterized by two index sets C, V satisfying the following condition: 
Condition I 

1) C,VcM, 

2) £np = 0, 

3) D/0, 

4) C = only if \V\ = 1. 

For simplicity, we will denote queue by Q^, and queue Q%\ by Qi- Also, we use the notation p!fc to denote 

a packet that is stored in virtual queue Qjy. The policies ensure that the following properties holds for any packet 

peQ%: 

Basic Properties of packets p G Q^: 

1) Any receiver % G V is a destination for p and any i G C is a listener of p. 

2) If p contains a native packet for user i T>, then this native packet has already been decoded by (i.e. is 
known to) i. Equivalently, p contains no unknown native packet for any user in V c , where c denotes set 
complement. 

3) For any native packet q for user i that has not yet been decoded by i, there exists exactly one packet p G Q£ 
(for some C, V) such that q = p(i). 

Hence, each packet p G Qffy can be written as p = ffij<=x>p(«) © d. We note that these queues are virtual in the sense 
that there are N real queues Qf at B, where queue Qf, i G N, contains native packets intended for user i that are 
unknown to this user. When a packet of flow i G Af arrives at B, it is placed in queue Qf. The virtual queues can 
be tracked by maintaining appropriate pointers to the real queues. Hence, at any instance in time, we can partition 
the real queue Qf into disjoint queues Q%(i), i G £>, where packet q G Qf is also in Q%(i) if for some p G Q^, 
it holds q = p{i). Therefore, viewing each queue as a set of packets, we have 

Qf = U QvW 

c,v 

where the union is over all sets C, V satifying Condition I. 

The virtual queues describe in effect the information that is gathered (through feedback) during the system 
operation for each of the packets in these queues, and are helpful in the description of the policies to be presented. 
The real queues Q^{i) will be useful in the analysis performed in Sections |Vl IVT1 We classify virtual queues into N 
levels, where level w G {1, . . . , ./V} contains all queues Q£ such that \C U V\ = w. Moreover, we classify queues 
of level w > 3 into sublevels, where sublevel w.u includes queues of level w with \C\ = u, u G {1, . . . ,w — 1}. 
In Figure \T\ we give an example of virtual and real queues when N = 3. 

In general, during the operation of the system, XOR combinations of packets are transmitted, at most one from 
each of the queues Q^. While the specific choice of packets depends on the received feedback and the specific 
algorithm that is employed, the following rule always holds. 

Basic Coding Rule 
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Figure 1. Virtual and real queues. 



A set V = jp^* , • • • ,Pj) m I of m = \V\ packets, one from each of the different virtual queues ^Qx>\ , ■ ■ ■ , QjT }> 
can be combined (by XORing) into a single coded packet only if 

V n C C r , Vr / n, n, r G {1, . . . , m} . (1) 

In other words, the destination set of each packet is a subset of the others packets' listener sets. 

Note that the Basic Coding Rule implies that V n C\V T = 0, for all r ^ n, n, r G {1, . . . , m}. Indeed, i G £> n 
implies, through (Q]), that i G £ r and, since according to Condition I it holds V r n £ r = 0, it follows that i ^ XV 
In fact, Condition I and the Basic Coding Rule imply that any user i G U™ =1 P n that receives the transmitted XOR 
combination will be able to decode its native packet. 

Lemma 3. Assume that Algorithmic Features [7J [2] are satisfied at the beginning of slot t. Then, if the packet to be 
transmitted at slot t is composed using the Basic Coding Rule, Algorithmic Features [7J [2] are also satisfied at slot 
t+1. 

Proof: The first part of the Algorithmic Features (i.e. any transmitted packet is a XOR combination of native 
packets) is trivial so we concentrate on the second part and assume that the transmitted packet p = @T =1 Pv k 
contains some unknown native packet q for user i G M (i.e. q is contained in some of the packets p v k k )- By the 
Basic Properties, it must hold i G U^L 1 I > fe, otherwise q would be known to i. Since are also disjoint sets, 
it follows that there exists exactly one k* G {1, . . . ,m} such that i G and q is contained in . But now 

i G T>k' is a destination for pJf t = q © c, where q is the unknown native packet for i and i is listener of c. The 

k * r r 

Basic Coding Rule implies that i G C r , for all r ^ k*, so that it holds p = g©c©@ r , k , p v r . Finally, p v k * t cannot 

c k * 

contain an unknown native packet for i other than q, since i is a destination for p^ k ,- Hence, i is a destination for 
p and the proof is complete. ■ 

Next, we describe how packets are moved between virtual queues (and, consequently, between the corre- 
sponding real queues Q%(i)), based on the received feedback. Generally, packets may be placed in various queues 
such that the Basic Properties are satisfied. For example, assume that packet p = p\ 2 © p\ 2 is transmitted (this 
combination satisfies the Basic Coding Rule), and that only user 2 receives the packet; hence packet pf 2 (2) reaches 
its destination and is removed from the real queue Q\ 2 (2). For the other packet movements, two choices are 
consistent with Basic Properties: 

1) Packet p\ 2 is moved to virtual queue Qf" and packet p^ 2 is not moved - hence, regarding the real queues, 
only packet p\ 2 {^) lS moved to queue Q 23 (l). This is consistent with Basic Properties since, after receiving 
p, receiver 2 becomes a listener for p\ 2 = p © p^ 2 , while receiver 3 is already a listener for that packet. 
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2) Packet p is moved to queue Qf 3 and packet p\ 2 , P3 2 are removed from queues Q\ 2 , Q3 2 respectively - hence, 
packet Pi 2 (l) is moved to queue Q? 3 (l) and p 3 2 (3) is moved to queue Qf 3 (3). 
Intuitively, the higher the level (and, within the same level, the higher the sublevel) of a virtual queue in which a 
packet p is stored, the better are the chances of combining the packet with other packets in such a way that multiple 
successful receptions are effected with a single transmission. Hence, when multiple choices for packet movement 
arise, we select the one that ensures that all packets involved in a transmission are placed in a higher level (and 
within the same level higher sublevels), else they are not moved. Hence, in the example above, we choose the first 

option. This is the underlying rationale on which packet movement is based in the following description. 

r r 

Rules for Packet Movement: Let packet p of the form p = p^ © . . . ®p^ satisfying the Basic Coding Rule be 
chosen for transmission, and let 5 be the set of users that receive p (the packet is erased at the rest of the receivers). 
If 5 = 0, the packet is retransmitted. Otherwise, 5/0 and let 5 be the subset of receivers in 5 n (U^L 1 £fc) 
with the following property: i G 5 iff it belongs to at least m — 1 of the sets C k , k G {1, . . . ,rn}. Hence, before 
transmission of p, user i G 5 is a listener for all but at most one of the packets pjf k G {1, . . . , m}. 

Every receiver in 5 which is a destination of p can decode its corresponding packet i.e., for any k G {1, . . . , m} 
and i G T> k D 5, packet p£ fc (i) is delivered to destination i. This is performed as follows: the Basic Coding Rule 

^ CP 

implies that i G C r for all r / k. Hence, this i G T> k H 5 can compute p v k k = p © ® r7 ^ Pjf , where all RHS 
terms are known to i after transmission of p (p is received correctly by i, since i G 5). But i G T> k n 5 is also a 
destination for p^J" , so that it can decode the unique unknown native packet p^ k k (i) contained in p^ . Notice that 
any i G T> k n 5 becomes a listener for p^ k after receiving p. 

If V k — 5 = 0, i.e. all users in V k receive p, then packet p v k is removed from virtual queue Q v k . Hence, if 
U'l n = l'D k — 5 = all packets are delivered to their destination and are removed from the corresponding queues. 
Otherwise, \J^ = ^D k — S ^ and we distinguish the following cases: 

1) it holds S = S - U^ =1 (£ fe U V k ) = 0, equivalently, S C U^ =1 (£ fc U V k ). Notice that, for m > 1, the 
latter condition is equivalent, by the Basic Coding Rule, to S C U^Li^, while for m = 1 it reduces to 
S C C\ U In both cases, and for each k G {1, . . . , m}, packet p v k k , where Vk — S 7^ 0, is moved to 

queue Q^ k ^g knS,)US and each packet p^{i), i G T) k — S, to queue Q^ L ^5 >fcn<S ^ U ' S (i). The consistency of 
these packet movements with Basic Properties is subsequently proved in Lemma [4] Hence, according to this 
rule, a packet is either not moved (if V k n S = 0), or is moved to a higher level (or within the same level 
but higher sublevel) queue. 

2) it holds S = S- U^ =1 (4UD{;) ^ 0. Again, this condition is equivalent to S = S- W% =1 C k / 0, for m > 1, 
and S = S — (C\ U T>\) ^ for m = 1. We further distinguish two subcases: 

a) If |(n£L x £fc U S) U (Uj™ =1 P fe - 5)| > nux fc= i,..., m \C k U P fc |Qthen packet p is moved to QC!^-5' 
packets are removed from queues Q v h and for each i in V k — S packet p^>* (i) is moved from 

queue Q^ k {i) to queue Qym" 1 ^^ (i). Lemma |4] shows again that this packet movement is consistent 
with Basic Properties and the packets are moved only to higher level or sublevel queues. 

b) If |(n™ =1 £ fc U S) U (U^ =1 P fc - 5)| < max fc=1 m \C k U V k \ then 

i) if 5 n (U^ =1 (£ fc U V k )) = , no further action is taken. 

ii) else, set 5 <— 5 n (U^L 1 (£fc U T> k )) and follow the rules of case[T] 
The validity of the above actions is proved in the following result. 

Lemma 4. Assuming that the Basic Properties are satisfied at the beginning of slot t, then the application of the 
Basic Coding Rule and Rules for Packet Movement to the packet transmitted at slot t satisfies the Basic Properties 
at the beginning of slot t+1. 

Proof: Let the transmitted packet p at slot t have the form p = ^ k=1 p^ k k where the various C k , T> k satisfy 
the Basic Coding Rule. We use the notation of the Rules for Packet Movement and examine each case of the 
Rules separately. The case \J^ = ^D k — 5 = (equivalently, \J£ = {D k C 5) is trivial, since no packets are moved 
between queues (hence, the Basic Properties still hold at t + 1) and all users in \J^ = ^D k decode their unknown 
native packets, as explained in the Rules for Packet Movement. We now assume \J^ = ^D k — 5 7^ and, for brevity, 

'it is easy to verify that this inequality is alwaystrue for m = 1. 
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only examine in detail the case m > 1, since m = 1 can be handled as a special case. Hence, we assume m > 1 
and further distinguish cases as follows: 

1) It holds S = S- U^ =1 (C k UV k ) = 0, i.e. S C U^ =1 (4UP fe ) = Uj^ =1 £ fc . In this case, for any k = 1, . . . , m 

such that T>k — S / 0, the packet is moved to virtual queue Q^ k ^\? knS ^ uS so that we need to show the 
following properties (shown in bold face) to be true at the beginning of slot t + 1: 

• All users in the non-empty set T>k — S are destinations for p£ fc and all users in U (V k fl5)u5 
are listeners for p^*: the destination property follows immediately by the Basic Properties at slot t, 
since T> k — S C For the listener property, we examine each set in the union Cu U n S") U 5 
separately: again, by the Basic Properties at t, all users in are listeners of p^\ Similarly, any user 

i G Pfc l~l 5 knows the value of p^ as well as the native packet intended for i (see the explanation in 

k r 

the Rules for Packet Movement). Notice that p£, fc cannot contain a second unknown native packet for 
i, since i G V k n S is also a destination for p^ k . Hence, any i G V k n 5 is a listener of p^ fe . Finally, 
consider any user in 5 n fl H <S) C = S n fl X>£ (since it holds S Q S) where we concentrate on 
users which do not belong to C k U (Vj, n <S). By construction of S, any user i G <S fl C k n belongs to 
all sets £ r , for r 7^ At, since it does not belong to and it must belong to at least m — 1 of the C sets. 
Hence, since this user received p, it can deduce p^£ = p © r ^fc • Furthermore, since i G P£, the 
Basic Properties at t imply that p^£ cannot contain an unknown native packet for i, so that i is again a 
listener of ptjf . 

• Packet p^ k k contains no unknown native packets for any user in set (V^ — S) c : consider any user 
i G (Vk — S) c = V c k U S = T> c k U (S fl T>k). The desired property now follows immediately from the 
Basic Properties at t (i.e. pjf contains no unknown native packets for any user in V k ) and the fact that 
all users in n S become listeners for p^ k k after successfull reception of p. 

• For each native packet q that has not been decoded by its intended user, there exists exactly one 
packet p G Qj), for some £, V, such that q = p(i): to prove this property, we use the fact that it is true 
at the beginning of slot t and, after transmission at slot t occurs, the only native packets that are decoded 
at the end of slot t are those intended for users in U^L 1 Pfc fl S. Hence, it suffices to only consider any 
unknown native packets contained in p^ , for all A; G {1, ... , m} with T>k — S / at the beginning of 
slot t. These native packets remain unknown to their destinations at the beginning of slot t + 1; however, 

due to the movement of p^ k k (i) to Q^ ) k k U( ;^ >knS ^ uS (i), these native packets are still included in exactly one 
packet, namely p^ k ^g knS ^ uS , and the proof is complete. 

2) It holds S = S- U^ =1 (4UDfc) / 0, equivalently S = S- U^ =1 £ fc . We only examine case 2.(a) of the Rules 
of Packet Movement, since case 2.(b) reverts to case 1, which has already been dealt with. Specifically, we 
assume |(n^L 1 £fc U S) U (U^ =1 X , fc — <S)| > maxj = i v .. ]m \Ck U T>k\ so that the composite packet p is moved 
to Qym" 1 ^^. Note that, applying verbatim the argument used in Case 1, we can show that all users in 
U^! =1 Pfc n S decode their unknown native packets. We now need to prove the following for the beginning of 
slot t + 1. 

• All users in U r k n =1 'D k - S are destinations for p and all users in rC^^Ck U S are listeners for p. For 

the destination property, consider any i G U^! =1 Pfc — S. Then, there exists some k* G {1, . . . , m} such 
that i G T>k* — S and, by the Basic Coding Rule, i G C r for all r / At*. Since % G T> k * is a destination 
for pfj£\ (Basic Properties at t) and we can write p^** = q © c, where q is an unknown native packet 

for i and i is listener of c, we conclude that p = q © c © @ r / fc „ p^, r where all RHS terms (except q) 

c c r 

are known to i. Since neither p v k l , nor any of p v r can contain any other unknown packet for i (recall 

that i is a listener for all p v r with r / At*), it follows that i is a destination of p. To prove the listener 

property, consider any i G n^L 1 £fe U S. If i G n^L^fc, then i already knows the value of p since it 

knows each individual component p^ , for At = 1, . . . , m. Furthermore, no p^ can contain an unknown 

native packet for i since i is a listener for all of these packets (due to Basic Properties at t); hence, i is 

a listener of p. It remains to consider the case i G S — n^L 1 £fc (notice that any such i learns the value 

of p upon reception). For any i G S — Pi^^Ck, i.e. i G S fl (U^L^^,), there exists some At* such that 

i G 5 fl But the Basic Coding Rule implies that 3P r 4 C V C T for all r ^ At* so that 
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r 

i € S n C k , n X>£. By the Basic Properties at t, no packet p~ with r ^ k* can contain an unknown 
native packet for i, since i € Hence, any unknown native packet for i can only be contained in "p-rf , 
provided that i € V k * (if i T> k *, the Basic Properties at t again guarantee that cannot contain 
an unknown native packet for i). But then, the conditions i G P^.and i £ 5 imply that i has already 
decoded its native packet, so that we can eliminate this case and conclude that any % € S — Pi™ =1 £k is a 
listener of p. 

• Packet p contains no unknown native packets for any users outside set U^ =1 Dfc — S. Consider any 
i € (U^L 1 Pfc — S) c = S U V^ k n = {D k . The property now follows immediately from the Basic Properties at 
t and the fact that all users in S become listeners of p after receiving it. 

• For each native packet q that has not been decoded by its intended user, there exists exactly one 
packet p 6 Qffy, for some £, V, such that q = p(i): the proof is similar to Case 1 by considering the 
users in \J^ = {D k n S and examining the native packet movement from Q^ k k {i) to Qy™ =1 ^-sW- 

■ 

Example 5. Suppose packet p = p\ 346 © plf © P3 246 is transmitted, so m = 3 and £>i = {1} ,£> 2 = {2,4} ,T> 3 = 
{3} , A = {2, 3, 4, 6} , £ 2 = {1, 3, 5} , £3 = {1, 2, 4, 6}. 

• Suppose p is received by users 2, 5 and 6, so S = {2, 5, 6}. It holds S = S - U 3 k=1 {£ k U V k ) = {2, 5, 6} - 
{1,2,3,4,5,6} = 0, so we are in case [TJ We have «S n (u| =1 (AUP fc )) = 5 = {2,5,6} and 5 = {2,6} 
because user 5 does not belong to m — 1 = 2 sets but only to set £ 2 . The 3 packets are moved as follows: 

- packet j> 2346 is not moved because V^nS = {l}n{2, 5, 6} = (equivalently, it is moved to Q^ lL ^? in5 ^ u<5 , 

i e q{ 2 > 3 > 4 ' 6 } u0u { 2 > 6 } = Q2346) 

- packet 'plf is moved to Q^n S) us ■ ^ ^^({^Ae})^ = Q i 23 56_ 

- packet pg 246 is not moved because T> 3 nS = {3}n{2, 5, 6} = (equivalently, it is moved to Q^_g 3 US , 

i e Qi^M^VM = Q1246) 

• Suppose now that p is received by users 7 and 8, so S = {7,8}. It holds S = S - U 3 k=1 (£ k U V k ) = 
{7,8} - {1,2,3,4,5,6} = {7,8} / 0, so we are in case |2] We have 

\{r\l =l c k \js)\j{\jl =l v k -s)\ = 

|(({2,3,4,6}n{l,3,5}n{l,2,4,6})U{7,8})U(({l}U{2,4}U{3})-{7,8})| = |{1, 2, 3, 4, 7, 8} j = 6. 



We also have 



max |A U2?fc| = max {({1,2,3,4,6} I , |{1, 2, 3,4, 5}| , |{1, 2, 3, 4, 6} |} = 5. 

fc=l,...,3 



Therefore, we are in subcase (a), and p is moved to Qya -1 ^^, i.e. Qi 2 34- 

If p is received by user 7, then S = {7}. It holds S = S - U 3 k=1 (£ k U D fc ) = {7} - {1,2,3,4,5,6} = 
{7} ^ 0, so we are in case We have | (n| =1 £ fc U <S) U (u| =1 D fc - <S) | = |{1, 2, 3, 4, 7} | = 5 and 

max fc= i 3 I £ fe UV k \ = 5. We also have S n (u 3 k=1 (£ k U V k )) = {7} n {1,2,3,4,5,6} = 0, therefore 

we are in subcase (b.i) and no packets are moved. 

If pis received by users 2 and 7, then S = {2, 7}. We have S = S -U 3 k=1 (£ k UV k ) = {2, 7}-{l, 2, 3,4, 5, 6} = 
{7} ^ 0, so we are in case El We have \(ri 3 k=1 £ k U S) U (u 3 k=1 V k - S)\ = |{1, 2, 3, 4, 7} | = 5 and 
max k=lt ... >3 \£ k UV k \ = 5. We also have 5 n (u| =1 (A U V k )) = {2, 7} n {1,2,3,4, 5,6} = {2} + 0, 
therefore we are in subcase (b.ii). Next, we set S <— <Sn (u k=1 (£ k U T> k )), i.e. S ^— {2}, and follow the rules 
of caseQ] We have S = {2} and the 3 packets are moved as follows: 

- packet p 2346 is not moved because T>\ n S = {1} n {2} = (equivalently, it is moved to Q^,_^ inUS , 

i. e .Q{W,6}U0U { 2} = Q23 46 ) _ 

- packet p\f is moved to Q^ nS ^ , i.e. J 5 ^« 2 ' 4}n{2 » U{2} = Q? 35 . 

- packet pg 246 is not moved because D3 n S = {3} n {2} = (equivalently, it is moved to Q^ 3 _J' 3 S , 
i.e. Qgf 4 .6}u0u{2} = g i 246) _ 
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As mentioned above, the rules for packet movement specified above are not unique. Also, information (who knows 
what about each packet) may be lost with this movement (see case 2.(b)). 

As explained, this choice is made on intuitive grounds in order to keep the system manageable and amenable to 
analysis. However, as will be seen in the next Section, for N = 4 even a more restrictive choice of rules suffices 
to implement a policy with asymptotically (as packet length increases) maximal stability region when the channel 
erasure probabilities are i.i.d. 

We next present two Lemmas, which follow from the Basic Coding Rule and Rules for Packet Movement, that 
will be useful in Section IVIII Lemma [6] states that a packet p, which is an XOR combination of packets that are 
stored in queues of at most level k, can include at most k such packets. 

Lemma 6. Consider packet p = pfj£ © ... © p^ m , where \T>i U £j| < k, i € {1, . . . , m}. Then m < k. 



m m 



Proof: It holds \V\ U£ a | < k. Also, since D T C C\, Vr G {2, . . . ,m}, we have (J V r C C\ and [j V r C 

r=2 r=l 

Vi U C\. Since n Vj = 0, i ^ j, it holds 



in 



< |X>i U jCi| < fc (2) 



Since V r ^ 0, r G {1, . . . ,m} , i.e, \V r \ > l,we conclude that J2T=i l-^Vl > m > which together with © implies 

m < k, as desired. ■ 

r r 

Let packet p = p^ © . . . (Bp v m be transmitted. Then, depending on the feedback from the users and in accordance 

1 m c c 

to the Rules for Packet Movement, either p is moved as a whole to a virtual queue, or packets p^ 1 , . . . , p v m are 

moved separately to virtual queues. More precisely, the following Lemma follows immediately from the Rules for 

Packet Movement. 

Lemma 7. After transmission of a packet at slot t, let packet p (not necessarily the transmitted packet) be placed 
at a virtual queue of level k.n. Then, either a) p is a combination of packets that at the beginning of slot t were 
at queues of level less that k, or b) p is a copy of a packet that at the beginning of slot t was either at level r, 
r € {0, . . . , k — 1}, or at sublevel k.l, 1 < I < n — 1. 



IV. Stabilizing Scheduling Policy 

In this Section, we investigate the design of policies that, under the coding restrictions and packet movements 
described in Section [TTTl stabilize the system whenever possible. We first need some definitions. 



A. System Stability and Stability Region 

Let X (t) , t = 0, 1, . . . be a stochastic process. 

Definition 8. [Stability] The process X (t) , t = 0, 1, . . . is stable iff 

lim lim sup Pr (X (t) > q) = 0. 

Consider next a time-slotted system U. At the beginning of each slot, a number of new packets belonging to a 
set Af of "flows" arrive to the system. Newly arriving packets of flow i € M are placed at infinite size queues, 
i.e. no incoming packets are ever dropped. These packets are processed by a policy it belonging to a set II of 
admissible policies. During this processing, packets are moved between a set of queues Q, and, at some time, 
they exit the system. Under policy it, let Qf (t) be the number of packets in queue Qi G Q at time t, and define 
Q n (t) = 5^Q,eQ Qf (*) • L et a l so Ai (t) i € jV be the number of flow i packets arriving at the beginning of slot t. 
For the purposes of this paper, we assume that the process {A(t)}^ , where A (t) = {A (t)} i£ j^f, consists of i.i.d 
vectors with A = E [A (t)\ . 

Definition 9. [System Stability] 

1) For a given arrival rate vector A, system hi is stable under policy it if the process Q^{t) is stable. 
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(R,x,x) 




Figure 2. Possible movements of packets p?2(l), P?2(2), p-j 2 (3). Destination of user i is denoted as di. Received feedback is denoted as 
(ui,U2) ua), where m is the feedback from user i, where R, E stand for received, erased, respectively, while X denotes an indifferent value 
(either R or E). 

2) The stability region of a policy, 7t n , it £ II, is the closure of the set of arrival rates for which hi is stable 
under it. 

3) The stability region TZjj of system hi under the set of policies IT is the closure of the set U^gn^- 

4) A policy vr* € II is stabilizing if TZ U = TV" . 

Consider now the system under study in the current work. At the beginning of each slot, a decision must be made 
at the base station concerning the combination of packets from the virtual queues that must be XORed to form the 
packet p = p v ^ © . . . © p v m to be transmitted. Such a decision is called a "control" -^''"'-p™ and we denote the 
set of such controls by I. Notice that, by definition, a control is identified by the set of pairs {(£>«, A)}™, and not 
by their order, i.e., control I^'"'^ is identical to control 1^"^'''"'^^' for any permutation a(i) of the indicies. 

We assume henceforth that the Basic Coding Rule is followed for the formation of packet p. For this system, an 
admissible policy consists of selecting, at the beginning of each time slot, one of the available controls '"' v m to 
form a packet p for transmission. After p is transmitted, packets are moved among the real queues Q^ k (i) according 
to the Rules for Packet Movement described in Section JII] 

A characteristic of such movements is that the destination of a packet movement is random and determined by 

3 12 

the feedback received after packet transmission. For example, assume that N = 3 and control I 12 3 is taken, i.e. 
packet p = p\ 2 ©P3 2 is transmitted. The packets involved in this transmission are pf 2 (l), ^12(2), P;P(3). Figure [2] 
shows the possible movements of these packets according to the received feedback. 

Our objective in this Section is to design a stabilizing policy for the system described above. It turns out that this 
system falls in the class of systems whose stability has been studied in 12. We next summarize the formulation 
and main results in [1] that will be useful in the development that follows. 

Consider a slotted-time network with a node set M. U d and edge (i.e. link) set £, where the special node d 
represents the destination of traffic originated at the nodes in M.. The node set M. is identified with the set of 
all real queues Q%{k), with k G V, i.e. each node m € M. corresponds to exactly one real queue. Let £^ t , 
denote, respectively, the set of outgoing links and incoming links to node m. We allow self-loops in the network, 
i.e. for node m € M., there may be a link (m, m), implying that the sets S^, both contain node m. A finite 
set of controls X is available. For each control / € X, "transmission" takes place over the sets of outgoing links of 
node m£jV, £™ t , as follows. 

• If, at a given slot, control I £ X is applied, then for any node ra G M at most /t m (I) > packets may be 
transmitted over the set £™ t in the following random manner: 
- For each m € M. and I £ X, there is a random sequence {R™ where each i?™ (I) takes values 

in the set £^ t , with the following interpetation. The n-th packet transmitted over the set £™ t when control 
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/ is applied, is received only by the recipient of the link R™ (I) . Of course, if R™ (I) = (m, m) then 
the packet is not received by any node in , hence it remains at node m. 

For a given n and /, the random variables {R l n (I)}.^ may be arbitrarily correlated. Moreover, we assume 
that for each control I £ X, the random sequences {{-R™ (/)} me ./vj} _ x are i.i.d and we denote p™ (I) = 

Pr (iC (/) = e) , e G £™ t (I), so that 

E *ro = 1 (3) 

To describe the stability region Tin of this network, we need some preliminary definitions. For control I £ 1, 
we define the set T(I) of vectors /, 

T (/) = {/ = {f e } ee£ : fe = P?(I)Vm, < /x m < £ m (I), » € AT, e G £™ } , (4) 
and the convex hull 7{ of the sets r(J) as 

ft = conv (r (I) , I G X) (5) 

The stability region of the network is described by the following Theorem. 

Theorem 10. [I] The stability region TZu of the system is the set of arrival rates A = {^m} m cu, A m > 0, for 
which there exists a vector f G H such that for any node m G A4 it holds, 

E fe + \ m < E /«• ( 6 ) 

We will apply the formulation described above to the network consisting of the real queues Q^(i), i G P, i.e., 
we consider A4 = (Q^(i)i * £ F° r tms network, since only one packet is transmitted per slot from any queue 
m, we have p, m (I) = 1, m G At. Also, the packet transition probabilities p™ (I) can be easily calculated. An 
example is given below. 

Example 11. Consider the case N = 3 and assume that control ij^ 3 is chosen, hence a combination p = pf 2 ©P^ 2 
is transmitted, where pf 2 = P12 (X) ® P12 (2) an d P3 2 = P3 2 (3)- The transition probabilities are then as follows: 

. Packet p\ 2 (1). 

1) If p is received by user 1, p| 2 (1) i s removed from Qf 2 (1) an d delivered to di (to d for the equivalent 
network). This event has probability ^0,(1} • 

2) If p is erased at user 1 and received by user 2, packet p\ 2 is moved to virtual queue and p\ 2 (1) is moved 
to queue Qf 3 (1). This event has probability P{\} t {2}- 

3) If p is erased at users 1 and 2, p\ 2 (l)remains at Q12 (1) • This event has probability Pup},®' 

• Packet pf 2 (2) . The transition probabilities are determined as in the previous case, by interchanging the indices 
1,2. 

. Packet pf (3). 

1) If p is received by user 3, pf (3) is removed from Qf (3) and delivered to d%. This event has probability 

^0,{3}- 

2) If p is erased at 3, pf (3) remains at Qf (3). This event has probability -P{3},0. 

The only difference between the network (At U {d} , £) and our model is that, in the latter, there are N packet 
destinations, di, i G A/ (one for each of the receivers) instead of a single one. However, we can combine all these 
destinations to a single destination d, so that any packet arriving in dj is considered to arrive at d. This affects 
neither the admissible policies, nor the queue sizes at the various native queues at the base station. Hence, system 
stability is not affected, provided that we are interested in the total queue size at the base station. 

Two implementation issues are worth mentioning at this point. First, there must be a mechanism for the receivers 
to know the constituents of the XOR combination of each received packet, in order to be able to use this packet 
in the decoding process. The simplest way to implement this is to use packet addresses to identify the packets 
involved in the XOR combination of the transmitted packet. These addresses can be placed in the packet header. 
Reserving bits to describe packet addresses implies some loss of throughput due to the introduced overhead. To 
simplify the description, in the current and next Section we do not take the overhead into account and address 
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the issue of stability in packets per slot. In Section IVII1 we discuss the number of addressed needed and loss of 
throughput due to overhead. 

The second issue is that under the schemes described in Section JIIJ the receivers need to save received packets so 
that they can correctly decode at a later time. Hence, if we are interested in taking these queues into consideration, 
we must ensure that the system remains stable even if the sizes of these queues are added to the total queue size 
at the base station. In fact, if the receivers are never informed by the base station as to which of their received 
packets will not be needed in the future, it is easy to devise scenarios where the queue sizes at the receivers grow 
to infinity even though the queues at the base station are stable (notice that system stability is only defined with 
respect to transmitter side queues). A simple way to deal with this problem is described in Section IVIIj 



B. Stabilizing Policy 

Applying directly the results in 12, we obtain the stabilizing policy described below. At the beginning of each 
time slot, the policy chooses a control of the form I = I-p 1 '"''-^ € X , where all real queues in the set M. (I) = 

U (J (k) > are nonempty, and forms the appropriate packet to be transmitted on that slot, p = ©j™ . 

If control / is chosen, one packet from each of the real queues in the set A4 (I) may be moved to another queue 
inside the network, or may reach the destination (and thus exit the network). No packets from any of the other 
queues are moved. The algorithm for choosing the appropriate control is the following. 
Algorithm 1 At each decision slot : 

1) For each control I = Zp 1 ''"'^"* G X : 

• Form the weights 

d (I) = max < Ki - £ pi (I) Kj, I , i € M (I) , 

e={i,j}G£* ut J 

where Ki is the length of the queue corresponding to node i. 

• Form the reward under the given control, 

C(i)= £ a (I). 

ieM(l) 

2) Select the control that maximizes the reward, i.e. I* = argmax/ e xC (I) . 

V. Outer Bound on the Stability Region 

In this Section we derive an outer bound on the stability region of the system under study by deparameterizing 
(i.e. eliminating the flow variables / in) Theorem [lOl This bound is identical with the bound on the information- 
theoretic capacity region of the broadcast erasure channel with feedback presented in 0, (3)- Although it was 
shown in |8] that the capacity region of the system under consideration is the same as the stability region of the 
system, we cannot directly invoke this result to derive the stability region outer bound via the capacity outer bound 
in 121, 0. The reason is that the latter capacity bound does not take into account the case of slots without any 
packet transmission, i.e. idle slots, so that, in principle, coding algorithms may take advantage of idle slots to 
increase capacity beyond the outer bound in (H, (3j. To distinguish the two channels, we call the broadcast channel 
studied in 0, the "standard" broadcast channel, and refer to the channel under study in this paper (i.e. the one 
containing idle slots) as the "extended" broadcast channel. 

As will be seen, the capacity of the standard broadcast channel, measured in information bits per transmitted 
symbol, differs from the capacity of the extended broadcast channel by at most 1 bit; in fact, this difference 
decreases exponentially w.r.t. the packet length L. Specifically, the following Theorem is proved in the Appendix. 

Theorem 12. A capacity outer bound C ou t, measured in bits per transmitted symbol, for the N-user "extended" 
broadcast erasure channel with feedback is given by 



C out = { R ■ max 
o-eP 



El 2~ l ^aA<l\ (7) 



13 



1 



where V is the set of all permutations on Af and A a = YlkeAf l _ E{ — — . 

We first describe the stability region of Theorem \10\ in a form that is more convenient for calculations. Any / 
in H can be written in the form 

f = J2<t>if (I) , for some > 0, < 1, (8) 

lex lex 

where 

/(I) = (/e(/)) e6£ , 
f e (I) = p™(I)Vm(I), < fi m (I) < m G M, e G = 



and, for any control / = Ip 1 '*"^" 1 , 



Hence it holds, 



and 



Am GO - <! oth e r ^ S g } • ( 9 ) 



E /e= E E^w 

ee£™ t e6£™ t /eX 

= E^ E pT(i)^ m (i) do) 
£/.= EE^w 

ee£™ ee£> /ex 

= E E HViWpKI)- (ID 

/ex e=(J,m)e££ 



Replacing (fTOb . (TTTb in ©, we have 



+ A m <X)( 1 -P^m)( J ))MmW0i > (13) 
/ex 



Yl & (I)p l eW + A m < E PeWV™ (I) J , ™ G M (12) 

/eX \e=(Z,m)££!£ / /ex \ee£™ t 

or equivalently, taking into account ©, 

/ \ 

E mWdSJ) 

/ex I e=(Z,m)e£™ , 
\ l^m / 

Hence, the stability region 7£n of the system is described by either one of (fT2l . (TT3T ). combined with 

0<A*mGO<£mOO (I 4 ) 
</>/ > (15) 

E^/<! (16) 
/ex 

where fi m {I) is given by ©. Since, in our case, new packet arrivals are only placed in queues Qi{i) i G Af, it 
holds 

A = I ^ if ™> = Qi(i), i eftf 
m \ otherwise 

The next Theorem, which is proved in Appendix |A] describes the main result of this Section. 
Theorem 13. The following relation holds 

( 4 > 

-V(i) 



Tin ^ < A: maxV - — — — < 1 > = C„, (18) 



where V is the set of permutations on Af and a(%) = {a (1) , . . . , <r(i)}- 
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Table I 

Permitted controls for levels 1 to 4. 





Level 1 


Level 2 


Level 3 


Level 4 


Permitted controls 


Control 


Control 


Control 


Control 


h 


11 


Tjk,i 
1 i ik 


j]kl,i 
*i,jkl 




n 


P 
1 jk 


P 

1 jkl 






jjk,ik,ij 
i,j,k 


jkl, i] 
ij,kl 






jjk,ik 
i i 




rkl,tjl,ijk 
*ij,k,l 






pk 

i 


jkl 








jj kl,ikl,ijl,ijk 
i,j,k,l 








jjkl,ikl,ijl 
i,j,k 








j]kl,ikl 
i,j 








j jkl 

i 



VI. The Case of i.i.d. Channels: Stability Region for 4 Users 

In this Section, we assume that the erasure events for all receivers are i.i.d, and denote by e the probability of 
such an event. Hence, 

We consider the case of a channel with 4 receivers and show that, for all < e < 1, if A G C u , then A G IZn, 
i.e. TZu 5 Cu- Hence, in this case we have TZu = C u and the stability region using only XOR operations coincides 
(barring addressing overhead) with the capacity region of the standard broadcast channel, and is within one bit 
from the stability region of the extended broadcast channel (which, due to (§ ], is equal to the capacity region of 
the extended broadcast channel in Theorem [T2l) . 

To proceed, we restrict the set of available controls by allowing only intra-level coding, i.e. we only combine 
packets from queues of the same level. This restriction simplifies the calcuations and shows that even a restricted set 
of controls suffices to achieve the maximal stability region when channel erasure events are i.i.d. We note however, 
that if channel statistics are non-i.i.d., the additional controls are helpful in increasing the stability region. 

The set of permitted controls is described in Table Jl where i, j, k, I G {1, 2, 3, 4} are distinct. 

For the rest of this Section, we assume without loss of generality that 

Ai > A 2 > A 3 > A 4) (19) 

which implies that 

C^^Et^^ 1 } (20) 

We will show that if A G C u , then A G TZn, which is equivalent to solving the following problem for any < e < 1. 

Problem: If Ai > A2 > A3 > A4 and Ylt=i 1=71 < 1> And parameters (pi satisfying (fT3T)-(fT6l). where M is the 
set of all queues (i) , i G V, and C,T> satisfy Condition I. 

In the following, we will describe the procedure according to which \i m (J) , </>j, m G M., I G X, are calculated. 
First, we set 

Hm (I) = Am(-f) , m G M, I G X, (21) 

ensuring that (fl~4l is satisfied. It remains to determine (pi, I G 1. Notice that, for any given value of e, (1211 ) 
transforms (fT3T ), (fl~5T ), (fT6l ) into a linear program (LP) w.r.t (pi, so that achievability of the rate A is reduced to LP 
feasibility (a similar LP-based approach is used to describe an achievable scheme for a 2 user MIMO setting over 
broadcast erasure channels in Q). However, since e takes a continuum of values, we cannot solve the resulting LP 
for each e but need to determine (pi analytically. 
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To simplify the notation somewhat, for control I 



..,V„ 



we denote 



An overview of the approach follows. We start from inequalities ( fT3l ) referring to queues at level 1, i.e. Qi (i), 
and determine all fa , ensuring that these inequalities are satisfied. In general, having determined 4>i for all controls 
that involve queues at level I, we consider the inequalities ([T3l referring to queues at level I + 1 and determine 4>i 
for all controls that involve queues at level I + 1, ensuring that these inequalities are satisfied. During this process, 
it is ensured that ([TBI) is satisfied. After all 4>j are computed, it is checked that (fT6l ) is also satisfied. 

We now proceed with the detailed description of the manner in which , I £ 1 are determined. We will use the 
following terminology in the decription. If, under an allowable control /, it is possible to have a packet movement 
from queue i to j, we say that there is a "flow from queue i to queue j" under control / and we name Puj) {I), 
the "probability of flow" from % to j under control /. We also say that there is "flow from queue i to queue j" if 
it is possible to have a packet movement from queue i to queue j under some of the allowable controls. 
Level 1 At this level, there are 4 queues (equivalently, nodes in M.) of the form Qi(i) , i £ {1, . . . , 4}. There 
are no incoming flows from other nodes to Qi (i), but there are new packet arrivals of rate A, at every 
Qi{i). The only control that may result in packets leaving Qi {i) is /j, so inequality ( [TBI becomes 



Level 2 



Xi < (1 - e 4 



j. To satisfy this inequality, we set 



fa = Xi/(l-e 4 ), V»€{1,.. 

At level 2, there are 12 queues of the form Q\ {%) , i, j € {1, . . 
each of these nodes is under control /j, with probability e 3 (l - 
under controls 2|'j and If, that result in packets leaving with probability 1 — e 3 . Hence, inequality (TT3T ) 
becomes 



,4}. (22) 

, 4} , i 7^ j. The only incoming flow to 
e), while there are two outgoing flows, 



Similarly, for node Q % - (j) we have 



£ 3 (1 



< (1 



< (1 



3 ' 



(23) 



(24) 



Since fa, <j, 
parameters, 



•j have already been determined by 

(ftl'p fai, <P)- Next, assuming i < j we set 



there are now 2 inequalities with 3 unknown 



0, if i < j. 



(25) 



Intuition for this is gained from inequality ( fT9l ). Specifically, observe first that through ([8]), we can 
interpret 4>j as the long term proportion of time that control / is applied. Also, it can be seen that 
as long as both queues Qj (i) , Q l - (j) are nonempty, the algorithm in Section HV-B| chooses to apply 
control PA instead of I] or II, since control Ij'j potentially results in two packet transitions during 



h3 



3' 



one time slot, while If , Jj can lead to only one packet transition. During the time interval that control 
l\' % - is applied, the same number of packets is expected to have left queues Q\ and Q* , as erasure 
probabilities are identical. When one of Qj, Q % - is empty, then either I 1 - or If, respectively, is applied. 

Since we assumed Ai > A2 > A3 > A4 and i < j, it is more likely that Q l j empties while J|'j is 
applied, and this event becomes very likely as queue sizes increase. Therefore, in the long run, one 
expects that control Ij is never applied and this is expressed by (T25T ). The intuition for setting certain 
to zero in the following description follows a similar reasoning. To satisfy (l24l . we set 



'■j 



! (l-6).A i /(l-6 3 ).(l-6 4 ), i<j. 



We substitute 



tt^'j in d23l and, to satisfy (|23l , we set 



^ = (Ai-A i )6 3 (l- e )/(l-6 3 )(l-e 4 ), i<j. 



Observe that (1231 ), ( f27T > can be combined to 

^ = [A,,-A J ] + 6 3 (l-6)/(l-6 3 ) (1-6 4 



(26) 



(27) 



(28) 
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Level 3 



where [x] = max{a;,0}. It is easy to check that the quantities in (|26l)-([28T) are nonnegative. 
At this level, there are 12 virtual queues of type Q k - (corresponding to real queues Q k - {%) and Q k - (j)) 
and 12 queues of type Qj ( corresponding to real queues (i)), where i, j, k 6 {1, . . . , 4} ,i ^ j ^ k. 
• Incoming flow to (respectively, to both Q\- (i) and Q\- (j)) occurs under control with 
probability e 3 (l — e). Outgoing flows from nodes of this form occur under controls 7^' l | and I k -, 
with probability 1 — e 3 . While for each of the queues Q k j (i) and Q k j (j) there is one inequality 
of the form (fT3l) . it turns out that these inequalities are identical. Hence, for both queues Q k - {%) 
and Q k , (j) the following inequality holds 



% < a - 



#2 + C 1 



We set 



'j 



0, so inequality ( 1131) becomes 



Next, to satisfy 



we set 



e 3 (i - e )-#;<(i 



ij,k ■ 



ij,k 



6 3 (l-6).<^/(l-6 3 ), 



(29) 



(30) 



where the second part of the inequality only depends on e and A, by substituting 0v. from (126b - It 



follows that 



kk,ij 
ij,k 



> 0. 



Possible incoming flows to Q^ k are due to controls hJi^J^k^l^i^fk^ikj an( ^ P ossl ble outgoing 



a i . . < Ti.lfc jjkAkAk T jk'ik T jk 

flows are due to controls ij^i^iju **■• ' 



; - ,I] K , where k G {1, . . . , 4} ,i ^ j ^ k. For Qf{z) 
^ + e 2 (l-e)- + 



inequality (fT3l becomes 



e 2 (l 



< (1 



+ 

,2\ 



tJ- 1 



+ 



k.i 

' *r'i,k 
ijk,ik,ij 



+ 



+ ■ 



ijk , ijk,ik 



jk,i 



+ 



+ 



jjk,ij 
\k 



+ 



0^ 



(31) 



For Q lk and inequality (fT3l takes the form of (f3TT >. with the appropriate exchange of indices. 
Specifically, for Q l - k and Q % k 3 , we have the following inequalities, respectively 



e) 2 

< f:i 



jk,ik,ij 
ij,k 



+ 6 2 (l-6) 



+ 



+ 



ik,j 



e 2 (l 



+ 



< 1 



rjk,ik,ij 
i,j,k 



+ e 2 (l 

k,ij 



+ 



+ 



+ 



ijkAk 



ik,ij 
'j,k 



t1 + 



+ 01 + 



Jjk,ij 
\k 



+ 



likAj 



+ 



ifcj 



+ 



(32) 



(33) 



All 4> parameters in the left parts of inequalities d3"l"T ), (1321 and (|33l have already been computed. 



Therefore, the unknown parameters at this point are 



k%3 Assuming, i < j < k we set 



rjk,ik,ij 
i,j,k 



ijk.ik 



<f>f = 

4>i- 



o. 

0, 

= 0, 
= 0, 



if i < j < k, 



and to satisfy (1331 ) we set 



dk,ij 
°Lk 



iik,ij 
b 3,k 



and 



(34) 



ijk,ik,ij 
i,0,k 



6 2 (1 



Kk 



+ 



til 



+ 6 2 (1 



+ 01 + 



i<j<k. 



fj^k 
ik,j 



+ 



fijk 
jk,i 



/(I 



(35) 
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(36) 



Next, we substitute d35t into d32l and set 

0g' ifc = [ e 2 (l- e ) 2 (^ + ^ + ^) + 

^ - *) • + ^ + C + - < 2 ; - C - <tf 

to ensure that d32l is satisfied. Finally, substituting d36l ) into (124b . we similarly obtain 

# = [^ 2 (l-^) 2 (^ + 4i + ^) + (37) 

We have to ensure that the quantities obtained in (|35T)- (|37T ) are nonnegative. The formulas are getting 
very convoluted at this point but they are easily calculated as functions of the erasure probabilities 
and the arrival rates using symbolic computation packages. Using such a package (we used Maple 
13.0), it is easy to see that the resulting expressions are nonnegative. 

Example 14. By computing 4>{ k / k we find that 4>\ k / k = e 2 (A 2 - A 3 ) _ e0 _fjJ.Xl'+e+i - From ® we 
have A2 — A3 > 0. Additionally, — e + 1 > 0, so the numerator is nonnegative. Also, the denominator 
is positive, because — e 6 — e 5 — e 4 + e 2 + e + 1 = (l + e + e 2 ) (l — e 4 ), where both parentheses are 



positive. Therefore, (jj k - lk > 



Level 4 At level 4, there are 4 virtual queues of the form Q\j k (which corresponds to real queues Q\ - k (i), 
Qijk (i)> Qijk (k))> 6 virtual queues of the form Q\ l - (which corresponds to real queues Q k j (i),Qij (J)) 
and 4 queues of the form Q? kl (corresponding to real queue Q^ kl {%)). 
• Incoming flows to queue Q\j k are due to controls I k j l k , Ij'jj.'l and if^^' 1 ^, with probability 
e 3 (l — e), while outgoing flows are due to controls I i fl k , and l\, k with probability 1 — e 3 . We set 



b\- k = so that inequality (TT3l becomes 



To satisfy (1381 ). we set 



l,ijk _ 3/-1 ,\ ( jk,ij 1 jj,ik , ,i,jk , jjk,ik,ij\ 1 /-, _ 3 



,,,, - ^ - «) w:i + C + *JS + / (' - ■ ^ 

Incoming flows to queue Qg are due to controls I k $, igf, J*J* ijg, ijg, /jfjg?'* 

with probability e 2 (l-e) 2 , and /g' ife , l|f , igg, i*$ with probability e 2 (l-e). Outgoing 

r ta,ij T kl,ijl,ijk . j-fcj 



flows are due to controls Z^f', ifji^' 1 ^ and Z-^ with probability 1 — e 2 . Therefore, inequality ([T3l 




'J 



so that (l40l becomes now 
Similarly, for queue Q^,, inequality (TT3l becomes 

^ 2 (i - ^) 2 OS + C + C + C + 4fl + 45 + «$!! + <B ,ife + ^ ,Jfc ) + 
?f + <f + tfj + <!) < (1 - ' 2 ) + ^S' 



(40) 



(41) 



(42) 



(43) 
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All </> parameters in the left parts of inequalities d40l ) and d43l) have already been computed. 

Therefore, there are two inequalities with 3 unknown parameters: 4\j'ki' $tj % kY^ k anc ^ ^li^Y^- 
Assume for the rest of this case that % < j < k < I. We then set 

^kl,ikl = (44) 



Now, there are only 2 unknown parameters, so to satisfy (1431 ) we set 

$2 = [e 2 a - OS + C + C + 4* + $5 + 45 + $2 + <#* + ^S?'* 

+e 2 (i - £ ) + <f + + O/ (i - , 2 ) 

(45) 

and to satisfy (l4Qb we set 

^fP = w - 'i 2 (<• + + + 45 + 4;] + oji 1 :, + $i + *Sf* + C 

*&* + + + 4^)]/ (i - e 2 ) - <g 

(46) 

For queues of the form Qf° l , incoming flows are due to controls of the form Jp 1 ''"'^" 1 , i G 
Pi, |PjU£j| < 3, j £ {1,...,3}, as well as controls of the form I^ 1 ''"'^" 1 , i G V\, \D\\ > 
2, |Dj U Cj\ = 4, j G {1, . . . , 3}. Outgoing flows are due to controls of the form I v u '"' v m , i G 
£>!, = 1, |Dj u Cj\ = 4, j G {1, . . . ,4}, with probability (1 - e). Therefore, inequality (fl3T> 
becomes 

j¥* j'.fcA j>M« / 





0, Va ^ i, a, b,c,d G {i , fe, /} 



e (i- e ) 2 E<^+ E^f+ E <@ + 

e(l-<0( E*P + E «$g + E ^fP ' 

\j,k^t j,k,l^i j,k,l^i 

(i-^)(Cf*+^+C" + 

. ijkLikLijk , ,jklAjLijk . ijkLikl . ± jkLijl , 

Similar inequalities can be formed for Q 1 ^ 1 , Q^ 1 and Q\^ k . In general, if i < j < k < I, we set 

^bcd,acd,abd = ^ y ( 0) b,c) / fc) , a,b,c,de {i,j,k,l} 

' '^,a C d = o, V (a, 6) 7^ (i, , a,b,c,d E {i, j,k,l} (48) 

Therefore, in the inequality concerning Qff the only unknown parameter is ]* We 
convert this inequality to equality and set 4^ j'kl '^'^ to the value satisfying it. Then, we proceed 
to the inequality concerning Q % j* 1 , where we substitute ^ k ^ k }^ l >v k an( j se t a value for the remaining 
unknown, (fP^j^ 1 ' 1 ^ 1 , such that it satisfies the respective equality. Then, in the inequality concerning 
Q l h l we substitute <^^^,«jMjfc an( j an( j sgt ^jkljM tQ ^ ya ^ ue sat i s fying the respective 

i t r-- 11 u * ijkLikLijLijk ijkLikLijl j ijkLikl . rjn\. j j.7'W ■ • -i 

equality. Finally, we substitute yifkl '^iik a j (147 D and set (p\ m a similar 
way. Again, it is easy to verify through symbolic computation that all computed quantities are 
nonnegative. 

Finally, to ensure that (fT6l ) is satisfied, we calculate the sum of all flows, and find 

lex i=i 

Since, by assumption, it holds Ei=i t^t < 1; we conclude that E/ex^/ — 1> as desired. Hence, we have proved 
the following result. 
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Theorem 15. For the case of 4 users, and for i.i.d erasure events, the stability region of the system is given by 

Moreover, the policy it* € II described in Section \IV-B\ usins the XOR controls described in Tablets stabilizing. The 
stability region coincides with the information theoretic capacity region of the standard broadcast erasure channel 
with feedback, and is within one bit (actually, 0(2~ L I A ) bits according to Theorem 172] where A* = m&x a A a ) 
from the capacity of the "extended" broadcast erasure channel with feedback. The latter is equal to the stability 
region of the system under any coding strategy. 

VII. Implementation Issues 

A. Packet overhead 

As mentioned in Section [TV] for the proposed network coding scheme to work, every user must know the identities 
of all native packets that constitute a composite (i.e. XOR combination) packet it receives. Having this information, 
a user is able to decode the native packet destined for it. A simple mechanism that can be used to provide users 
with this information is equipping every native packet with a Packet ID, which consists of the packet's destination 
and a sequence number. If a transmitted packet is composed of m packets, then it contains in its packet header the 
m packet IDs. 

To compute the overhead bits needed to implement the above mechanism, we need to find the maximum number 
of Packet IDs that may be included in a packet that is placed in a virtual queue of a certain level. This is expressed 
in Lemma [16] below. In the following, when we say that a packet comes from level k (or exits level k) we mean 
that it is an XOR combination of packets placed in queues of levels 1 to A; (with at least one packet being in a 
level k queue). 

Lemma 16. Under the coding scheme of Section [777] it holds 

a) Any packet placed in queues at sublevel k.n, n = 1, 2, . . . , k — 1, k > 2, contains at most (k — 1)! packet IDs. 

b) Any packet exiting level k > 2 contains at most k\ packet IDs. 

Proof: We use induction on k to prove the Lemma. For k = 2, the Lemma follows immediately from the 
Rules for Packet Movement in Section [Till We now assume that the Lemma holds for levels 2 up to k — 1 and 
show that it also holds for level k. We first prove part a) of the Lemma by induction on n. 

Part a): If a packet p is placed in a queue at the lowest sublevel of level k, i.e. k.l, then according to Lemma [7] 
p comes from levels I < k — 1. Hence, according to part b) of the inductive hypothesis, it contains at most (k — 1)! 
packet IDs, so that part a) holds for n = 1. Assume next that part a) holds for all packets p placed at any sublevel 
from k.l up to k.n with 2 < n < k — 1, i.e. assume that all packets p in sublevels from k.l up to k.n contain 
at most (k — 1)! packet IDs. We shall prove that any packet in sublevel k. (n + 1) also contains at most (k — 1)1 
packet IDs. According to Lemma [7J for a packet p at sublevel k. {n + 1), one of the following two cases holds. 

1. Packet p comes from level I, where 2 < I < k — 1. Then, according to part b) of the inductive hypothesis, p 
contains at most (k — 1)! packet IDs. 

2. Packet p was placed before the current slot transmission at a queue in a lower sublevel of the same level, i.e. 
a sublevel from k.l up to k.n. According to the inductive hypothesis on n, packets in these sublevels contain at 
most (k — 1)\ packet IDs. Since Lemma [7J states that packets from lower sublevels are merely copied to higher 
sublevels, it follows that the maximum number of packet IDs they contain remains the same, so packet p at sublevel 
k. (n + 1) will also contain at most [k — 1)! packet IDs. Therefore, packets at all sublevels k.n, n = 1, 2, . . . , k — 1, 
k > 2, contain at most (k — 1)! packet IDs. This completes the proof of part a) of the Lemma. 

r r 

To prove part b) of the Lemma, consider a packet p exiting level k. This packet is of the form p = © . . .®p^ m , 
where each p v r belongs to a queue of at most level k, hence the maximum number of packet IDs p may contain 
is the sum of the packet IDs contained in packets p^' , r E {1, . . . , m} , which is at most m(k — 1)!, due to part 
a). From Lemma [6] it holds m < k, therefore any packet exiting level k contains at most k(k — 1)! = k\ packet 
IDs. ■ 

Up to level 4, the maximum number of Packet IDs that may need to be included in a packet is 4! = 24. Assuming 
a packet ID of 20 bits and packet lenght of 12000 bits, the overhead is approximately 4%. It can be seen that 
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the maximum number of Packet IDs needed increases dramatically with the number of users N and it is very 
important to address this matter as N increases. As a future work, various suboptimal policies that reduce the 
necessary number of Packet IDs can be investigated. For example, the transmitter may choose not to send packet 
combinations if the resulting packet header exceeds a certain number of Packet IDs. Another policy towards this 
direction could involve coding of packets only until a certain level. Specifically, for N users, only the virtual queues 
until level I could be created, where I < N. In case a packet is received by more than I users, additional receivers 
would be ignored and the packet would be placed in one of the level I queues. 

B. Queue stability at the receivers 

As mentioned in Section [IVl another problem that may arise is possible instability of queues at the receivers, 
where all packets received by a certain user are stored. A simple way to avert this possibility is to take advantage 
of the fact that when the queue sizes at the based station become empty, all packets formed during previous 
transmissions are not needed at the receivers. Therefore, we can let the base station inform all receivers when its 
queues become emtpy, by, for example, leaving a slot empty after a series of transmissions taking place when the 
queues are nonempty. Under this modification, using standard results from regenerative theory, it can be shown that 
the system is stable if and only if the total queue size at the base station is stable. 

VIII. Conclusions 

In this work, we presented a network coding scheme for the broadcast erasure channel with N multiple unicast 
sessions based on the coding scheme we proposed in [4]. This scheme is characterized by low complexity, as only 
XOR operations are allowed. Also, instant decodability, i.e. the ability of any user that receives a coded packet 
to instantly decode its own native packet, is ensured. While in (4[ the possible transitions of packets between the 
virtual queues of the coding scheme were determined only for the case of 3 users, in this paper we give general 
rules for packet movement. 

Furthermore, we assumed random packet arrivals and presented a stabilizing policy based on this coding scheme. 
We then derived a capacity outer bound for the system under examination. For the case of 4 users and i.i.d. erasure 
events, we proved that any rate within the capacity outer bound is stable under the proposed policy, therefore the 
stability region of the system is identical to the capacity outer bound of the BEC channel with feedback. 

Finally, implementation issues were examined, such as the increase of packet overhead as the number of users 
increases, which is due to the number of packet addresses needed to completely describe a coded packet. The exact 
number of addresses needed in the general case of N users was found to be N\. Future work could be aimed 
towards the development of suboptimal variations of the proposed policy that will require a smaller number of 
packet addresses,this reducing packet overhead. 

Appendix 

A. Proof Of Theorem\l3\ 

We need some preliminary definitions. Define the sets M\ = and Mi = {1, 2, . . . , i — 1}, for i G M with i > 2, 
as well as 

Qi = {Qt> (») ■■ i G V, and C,V - {i} C Mi} 

X i = : ^ £ £>, and C,T> — {i} C^}. 

Notice that 2j P\lj = for i / j. This is due to the fact that the existence of a control / G Zj nlj would imply 
that i G Mj as well as j G Mi, which is impossible. We also define the set Q^ in the subnetwork consisting of 
queues (i.e. each node is a queue, as described in Section ITVl) as follows: 

Q¥ = [Q% (i); ie V, and V, C C M} U {d}, 

Denote with C ou t{Qi) the set of all outgoing links in the cut [Qi, — Qi], i.e. 

Com (Qi) = {e=(i,j)eS: i G Q u j G Qf - Qi} 
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and notice that, according to the Rules for Packet Movement, the set of all incoming links to the cut is empty, i.e., 

C in (Qi) = {e = (i, j) € £ : % G Sf - Qi, j G Qi} = (49) 
To prove Theorem [T3l it suffices to show that, under (fT2l ) and (fT4l-(fT6l). it holds 



V A,; 



E < 1, (50) 

since the same argument can be repeated verbatim for any permutation <r(i), i G M. Summing ([T2l over all m G Q« 
and using (fT71 ) yields 

E ^ E E W (/)PL(/) + A, < E E E PeWUm (I) <t>I, V* G A/" , 

or, rearranging the terms, 

A ^ E I E IE Pe(Ihm(I) ~ E W (I)p l e(I)) J 0/- (51) 

But d49l ) and the construction of C OU i(Qj), Ci n (Qi) imply 

E I E p?(ihm(i)- E wcop&oo 
= E wWpIw- E w(j)^(j)= E 

e=(i,m)eC (Q») e=(i,m)eC in (Qt) e=(Z,m)GC(Q0 

Under any control ^'p!, 1 "'xT £ it holds 

E W (/)pi(I) = l-e^_ M . (52) 

e=(«,m)GC (Qi) 

This follows from the fact that, when I-p'-p^' '"'xT i s applied, it holds /i; (/) = 1 for / = and pn (I) = for all 
other queues in Q^f , according to (©. Moreover, under '"'zT > a native packet for user i is transferred to one 
of the queues in Q^ — Qi iff the packet is not erased by all users in M — Mi ■ 
Hence, (IBTI ) yields through (l52l 

Aj < E ( x ~ W-aO 0/ 
and, summing over all i G M, we conclude that 

However, since 2j nZj = for all i 7^ j, it holds Y2ieJ\fYliei $1 — ^iei^ 1 — ^ anc ^ ^501 is proved. 



B. Proof Of Theorem [72] 

We first need to establish some notation and prove a few intermediate results. We consider the "extended" 
broadcast erasure channel (BEC), where the transmitter has the option of not transmitting in a given slot (as 
opposed to the "standard" BEC that appears in the literature). This is equivalent to considering that the transmitter 
sends in this slot a special (null) symbol, denoted as 0. Hence, in information theoretic terms, given a standard 
point-to-point BEC with an input alphabet of X and output alphabet of y = X U {*}, where * denotes an erasure, 
the extended point-to-point BEC has input alphabet X' = Xu{0} and output alphabet y' = X'U{*} = XU{*, 0}. 
Since we consider feedback, we assume that, if the transmitter sends symbol 0, all users send as feedback back 
to the transmitter. Hence, at slot I, each user can send feedback Z G {ACK, NACK, 0} to the transmitter, where 
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ACK (resp. NACK) denotes a successfull reception (resp. erasure) of a non-null symbol, while denotes a null 
symbol transmission (and reception). 

The N user version of the extended BEC follows from a simple "vectorization" procedure. Specifically, let 
M = {1, . . . , N} be the set of ./V users and denote with Wi the message for user i 6 M. The transmitted symbol 
at slot I is denoted as X(l) (with X(l) G X') and we also introduce the shortcut notation X 1 = (-X"(l), . . . , X(l)). 
Furthermore, let YJ(Z) G y' be the symbol received by user i at slot I, while Zi(l) G {ACK, NACK, 0} is 
the feedback sent by user % to the transmitter at slot I. We can also define an auxiliary random variable Zi(l) G 
{ACK, NACK} that is independent of X(l) and all previously generated random variables (up to slot I) so that 
it holds 

' Zi{l) if X(l) ^ 
if X {1) = 



Zi(l) 



Notice that, for any z ^ 0, the events {Zi(l) = z} and {Zi(l) = z,F(l) = 1} are identical. We now introduce the 
following "vectorized" entities 

W {1J] = (W^ . . . ,Wj) 

Yl = (Y i (l),...,Y i (l)) 
Y {1)j] (l) = m),...,^)), Y l [hj] = (Y {1)j] (l),...,Y [ltj] (l)) 
Z [hj] {l) = (Z 1 (l),...,Z j (l)), Z l {1J] = {Z [hj] {l),...,Z m {l)) 
Z [lj] (l) = (Z 1 (l),...,Z 3 {l)) 

and use the shortcut Y = Y^^, Y l = YL n ^ (with similar interpretation for Z, Z l ). 

The subsequent analysis closely follows the approach in ifTUl . with some necessary variations due to the fact 
that Z(l) are X{1) are not independent. The following Lemma can be proved by straightforward manipulations of 
information measures. 

Lemma 17. Let A,B,C,D be discrete random variables. The following identities hold. 

1) Conditioning can be added to either part of mutual information: 

I(A; B\C, D) = I(A, C; B\C, D) = I(A; B, C\C, D) = I(A, C; B, C\C, D) 

2) Let B be independent of the joint ensemble (C, D). It then holds I(A, B; C\D) = I (A; C\B, D). 

3) Let D be independent of the joint ensemble (A, B, C). It then holds I(A; B\C, D) = I(A; B\C). 

4) Conditioning can be augmented by redundant condition, i.e. if the event {B = b} implies {C = Cb}, it then 
holds H(A\B, D) = H(A\B,C,D). 

5) It holds I{A; B\C) = I{A; B\C, D) + I (A; D\C) - I (A; D\B, C). 

We now consider an arbitrary code £ for the extended BEC with feedback (see fl8j for a detailed description of 
encoding and decoding functions of (C) and denote ir(l) = Pr(X(Z) ^ 0) and F(l) = I[X(l) ^ 0]. The following 
results, whose proofs can be found, respectively, in sections |Cj O of the Appendix, will be used. 

Lemma 18. For any rate H = (R\, . . . , Rn) that is achievable under <£, and for any j G M, it holds 

j n 

n zZ R ^H [ h « 1 )) + ""(OKI -Hi,..J})I(W [1J y,X(l)\Y\-} i ,Z l -\F(l) = 1)] + o(n) 



where h(-) is Shannon's entropy function. 

Lemma 19. For any rate R, = (Ri, . . . , Rn) that is achievable under <£, and for any j G M, it holds 

j n 

Y^Rk > (1 -e{W+i})X;(l -K{l))I{W [lt3 y,X(l)\Y\~] +1] ,Z l -\F{l) = 1) 



n 

k=l 1=1 
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Applying Lemma [18] for j — 1 yields 

nE *= lflfc < o(n) + £ fe(?r(/)) + £(1 - 7r(0)/(^ [1) ,_ 1] ;X(Z)|yj-)_ Z<~\ F(l) = 1) (53) 

where the second line was produced by using the inequality 

(l-7r(0)/(W[ij-i];X(0^ 

iFlG^ijj^Qiy^ (54) 
= (1 - tt(O) [l(% J _ 1] ;X(Z)|F|- j 1 ] ,Z / - 1 ,F(/) = 1) +I(Y l j - 1 -X{l)\Y l -]_ iy Z l -\F(l) = 1) 
and applying Lemma [19] for j — 1, to the first term in the last line of (l54l . Hence, we arrive at 

^Rkf- )<o(n) + - l - f>(vr(/)) 



n 

k=l 



n 

+ ^(l- 7 r(Z))J(y;- 1 ;X(Z)|Yj-j_ 1] ,^- 1 , J F(Z) = l) 



(55) 



z=i 

We are now ready to prove Theorem [12] We only consider the identity permutation (i.e. a(i) = i), since all other 
permutations are handled similarly. Summing (|55T ) for j = 2, . . . , N, applying Lemma [18] for j = N and summing 
the results yields after some manipulations (which involve a change of order summation between j and k) 

N jj / N \ n n N 

n E 13^—: * E - E wo) + Ed - E ^ x^KU > = !) 

fe=l i J V J=l U'-Jl / J=1 ; =1 J=2 



(56) 



(57) 



+ ]T(1 - 7 r(Z))/(W [1)iV] ;X(Z)|y z - 1 , Z l ~\F{l) = 1) + o(n) 

For notational compactness, we hereafter denote A = YljLi i_ e | — _ - ^ a l so holds 

L > H(X(l)\F(l) = 1) = /(X(Z); y^, Z w |F(0 = 1) + # (X(/)|y'-\ Z^ 1 , F(Z) = 1) 

N 

= ^/(^- 1 ; x(z)|y;-)_ 1] ,z i - 1 ,F(z) = l) + H{X(l)\Y l -\Z l -\F{l) = 1) 

iV 

> ^/(y/- 1 ;x(z)|yJ-)_ 1] ,z i - 1 !J F(z) = 1) + F(W [1>JV] ;X(Z)|y i - 1 ,Z'- 1 ,F(Z) = 1) 

3=2 

where the second line is derived by applying the chain rule over j. Inserting (1571 ) into (1561 ) yields 

iV 

n V * < V [Ah(*(l)) + (1 - vr(Z))L] + o(n) (58) 

The RHS of ([58]) is separable in terms of tt(Z) and its maximum can be computed via standard derivative arguments. 
In fact, the maximum in the RHS of ([58] ) is achieved for tt(Z)) = 1+ Il/a for I = 1, . . . , n which yields 

N E 

nV < raAlog 2 (l + 2 L ' A ) + o(n) = nL + raAlog 2 (l + 2~ L I A ) + o(n) (59) 

^ 1 - e {W} 

Dividing by n, taking a limit as n — > oo and using the inequality ln(l + x) < x, for any x > 0, yields 

N R 

yKk < L + 2 -L/A A m 
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Repeating the above procedure for an arbitrary permutation a on M produces 

N R 

y Kk < L + 2 -L/A„ A(r 

where A a = Y2k=i ^ ~ an( * snice tne ^ ast inequality must be true for all permutations a, the proof is 

complete. 

C. Proof of Lemma !7<Sl 
Fano's inequality implies 

3 

nY J Rk = H(W lhj] ) = I(W {l ^Yl- v Z n ) + o(n) (61) 



fc=i 



with 



I{W {lS X{l)\Y\^ v Z^,F{l) = l)Pr(F(l) = 1) 

i-i 



(62) 



(63) 



i=i 

n 

. . n 

Applying the chain rule twice with different order yields 

i(w [1)i]; z(0,x(z)|yj^ 

= /(iy [lji; z(z)|i^^ 

and since Z(l) is independent of all previous random variables given X(J), d63l yields 

Z(i)|r|rj], = /(W [lii]; X(0|F|rj], - I{W m ;X(l)\Y\-] v Z l -\ Z(l)) (64) 

Furthermore, since knowledge of X{1) implies knowledge of F(l), it holds 

IiW^-X^Y^Z 1 - 1 ) = I^W^X^F^Y^Z 1 - 1 ) 
= I(W [lS F(l)\Y\-] ] ,Z l - 1 )+I(W [lS 
Combining (l64l . (I65T ) yields 

/(^ [lj] ;Z(/)|rJ-j p Z / - 1 ) = /(^ [lj] ;^ 

-/(^^(oiyjrj],^- 1 ^^)) 

Defining the set -Ztiji = {Z[i,j] '■ ( i • • • > we can compute 

/(iy [lj1 ;X(/)|r;-j ] ,Z / - 1 ,Z(/)) = J] I{W {l ^X{l)\Y\-] y Z l -\Z(l) = z)Pr(Z(0 = z 

zeZ[i,j] 

= E ^^^(Ol^fc)],^- 1 ,^^ = z,F(Z) = l)Pr(Z(/) = z)Pr(F(Z) = 1) 

2 £-2[ij] 

= E I(W [lS X(l)\Y\-] ] ,Z l -\F(l) = l)Pr(Z(l) = z)Pr(F(l) = l) 

2 £2[i,i] 



(65) 



(66) 



(67) 



I(W [hj] -X(l)\Y\-l,Z l -\F(l)) 
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where we exploited the independence of Z(l) from all variables up to slot Z and used the facts that F(l) = 
implies X(l) = and EzeZx - Pr(Z(Z) = z) = 1. 

To manipulate the last term in (l62l) . we define the set Zn = {Znji : Zn 7^ (0, . . . , 0), (*,...,*)}. In words, 
2hji is the set of feedback vectors in which at least one user in {1, ... ,j} successfully receives the transmitted 
symbol and sends back ACK. Notice that, for any z ^ the event {Z^ ^(l) = z} implies full knowledge of 

ynj|(Z). It now holds 

= £ [^(^lyj-jpZ'- 1 ^^]© =z) - ff(^ [1 , Jl |y{-j p z'-\y [M (z) J z [1)i] (0 = z)] Pr(Z tlj1 (/) = z) 

= £ [F^jiyJ-jpZ'- 1 ,^© = 1,Z [M (Z) = z)- HiW^lY^Z 1 - 1 ^^)^^ = l,Z m (l) = z) 

zeZ[x tj ] 

xPr(Z [M (Z) = z)Pr(F(Z) = l) 



£ [H(W m \Y[-j v Z l -\F(l) = 1) - iJC^lyj-^Z'- 1 ,^© = 1,X(Z))] Pr(Z tl)i] (Z) = z)Pr(F(Z) = 1) 
(l-6 {li ... j} )(l-vr(/))/(^ [lj] ;X(/)|y| 1 J p Z'- 1 ,F(0 = l) 



(68) 



In the transition from the third to the fourth line of (I68T ). we used the event identity {Zm ji(Z) = 2} = {Zry](Z) = 
z,F(l) = 1}, which is valid for any z € Znji, while in the transition from the fourth to the fifth line we used 
the facts that Zhj] is independent of all variables up to Z (including i ? (Z), X(l)) and knowledge of y^j^Z), 
Zhj](Z) = 2; implies knowledge of X(Z) for any z G <2rijl- 

Inserting (168T ). (I67T ). (l66l ) into (|6TT >. via (l62l ). and using item [5] in Lemma [T71 produces 



it 



J>* < o(n) + £ [WiJ];^(OI>t"j]> Z, ~ 1 ) + C 1 - 'ii i})0 " 7r(Z))7(W [1)i] ;X(Z)|yJ-) p Z i - 1 ,F(Z) = 1) 

fc=i 1=1 

n 

< o(n) + J2[H*(l)) + (l-e {w} )(l -n{l))I{W [x ^X{l)\Y\-] v Z l -\F{l) = 1) 



1=1 



(69) 



where we used the inequality I{W {l)j] - F(l)\Y l ~^ v Z 1 ' 1 ) < H[F(l)) = Zi(tt(Z)). 
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D. Proof of Lemma [79] 

Performing similar manipulations as in the proof of Lemma [T8j we can write 

j n 

n^R k = H{W [hj] ) > I(W [lS Y^ J+1] ,Z n ) = ^/(^ lj] ;r [lj+1] (0,^(0l^ 1 +1] ,^ 1 ^ 

k=l 1=1 



[1J+1]' 

z=l 

n 



= E E JCW[ij]; ^[u+iicoi^lrj+i]. ^(0 = *) p ^(o = *) 

n 

= E E WnjlJ^Hj+ilCOI^Irj+i]^'"'^^) = • z ' i? ( / ) = l)Pr(F(/) = l)Pr(Z(Z) = z) 

n 

E E [^(w- [1Jl |y|rj +ip z'-i^(o = z,F(o = i) 
-^Mi^rj+iii^.^tij+iico^co = *,no = i)] p^f(/) = i)Pr(z(o = z) 

n 

E E [^ [1j] i^) +1] ,^- i ,z(o = 2; ,f(o = i) 



(1 - E^ 1 " vr(/))/(^ [lj]; Z'- 1 , F(Z) = 1) 



i=i 



(70) 



where we used again the independence of Z(l) from all other variables. 
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